Mosaic length and finite interact ion- range effects in a one 
dimensional random energy model 

S. Franz/ G. Parisi,^ and F. Ricci-Tersenghi^ 

^ Umversite Pans-Sud, LPTMS, UMR8626, 
Bat. 100, 91405 Orsay cedex, France 
^Dipartimento di Fisica, INFM-CNR SMC and INFN sez. di Romal, 
Universitd di Roma "La Sapienza" , 
P.le A. Moro 2, Roma 00185, Italy 
(Dated: February 2, 2008) 

Abstract 

In this paper we study finite interaction range corrections to the mosaic picture of the glass 
transition as emerges from the study of the Kac hmit of large interaction range for disordered 
models. To this aim we consider point to set correlation functions, or overlaps, in a one dimensional 
random energy model as a function of the range of interaction. In the Kac limit, the mosaic length 
defines a sharp first order transition separating a high overlap phase from a low overlap one. 
Correspondingly we find that overlap curves as a function of the window size and different finite 
interaction ranges cross roughly at the mosaic lenght. Nonetheless we find very slow convergence 
to the Kac limit and we discuss why this could be a problem for measuring the mosaic lenght in 
realistic models. 

PACS numbers: 05.20.-y (Classical statistical mechanics), 64.70.Pf (Glass transitions), 75.10Nr (Spin-glass 
and other random models) 
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I. INTRODUCTION 



The paradigm of 'random first order transition' or one-step replica symmetry breaking 
(IRSB) theory, provides an elegant^ framework to conceptuahze the phenomenology of liquids 
approaching the glass transition [1]. Unfortunatly, this scenario is strongly based on mean- 
field models {2] and mean-field- like approximations to liquid theories and cannot be taken 
literally in the application to real system. The main node that has to be untied to establish 
the IRSB scenario as a convincing theory for real materials, is how mean-field theory should 
be adapted and modified to take into account the finite range of interactions. Though a 
fundamental theory of glassy systems in finite dimension is presently lacking, proposals have 
been made that modify minimally mean field scenario to take into account the finite interac- 
tion range. In ref. 2[ Kirkpatrick, Thiurmialai and Wolynes developed a phenomenological 
theory, known as 'mosaic picture', where it is postulated the existence of a coherence length, 
that grows on lowering the temperature. Below that length the system behaves essentially 
as a mean-field glass, while it would cross-over to liquid behavior at larger scales. It results a 
theory where relaxation is dominated by activated processes stemming from the competition 
between interface tension and a bulk configurational entropy. The mosaic picture has been 
recently revived and deeply clarified by Biroli and Bouchaud ^|, who showed that while 
usual (point-to-point) correlation functions are insensitive to the possible growth of the co- 
herence mosaic length, it is possible to define different "point-to-set" correlation functions, 
able to reveal the growth of the mosaic length. In turn, the mosaic length has been related 
to the relaxation time of ordinary, time dependent correlation functions [sl. These papers 
prompted on one side numerical simulations on kinetically constrained glasses jol and on 
realistic glassy models 7|, on the other to theoretical calculations for models on trees and 
under the Kac limit [8] . These last models are the natural starting point for understanding 
the mosaic picture, since their local properties are well described by mean-field theory j^. 
In js] the study of point to set correlation function has allowed to derive a detailed picture 
relating the relaxation in the Mode Coupling regime for T > to the one in the mosaic 
regime for T < T^. The calculation, supposedly exact, concern the behavior of disorderd 
glasses in the Kac limit. In order to understand its relevance for short range systems, it 
is necessary to study the properties of convergence to the Kac limit for finite interaction 
range. It has been found in ref. p] that in standard Lennard- Jones supercooled liquids, the 
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transition from high to small overlap as a function of the box size is much smoother than 
one would expect from the mosaic picture. This poses the question of what behavior one 
should expect when the range of interaction is not large. 

In this paper, we address this question in a minimalistic finite dimensional model display- 
ing IRSB behavior in the Kac limit. The model is a one dimensional version of the Random 
Energy model extensively studied in the context of stochastic models for reaction diffu- 



sion equations and evolving populations ll|]. This has two main advantages: on one hand 
the Kac limit can be studied directly by probabilistic arguments, without having to resort 
to replicas or cavity techniques, on the other the model for finite interaction range can be 
studied exactly by transfer matrices. 

A recent paper addresses the problem of finite range corrections to the mosaic picture in 
a related one dimensional XORSAT model 12| . That paper concerns the zero temperature 



limit, while we concentrate on finite temperature properties. 

The organization of the paper is the following: in section [TTl we define the model. Sec- 
tion [TTll is devoted to the definition of the point to set correlation we study. In section HVl we 
discuss theoretical approaches to the computation of this quantity. In section |V] we discuss 
the results of exact computations with transfer matrices. Finally we draw our conclusions. 

II. THE MODEL 

In order to compare the behavior of finite range interaction systems with mean-field 
theories we need a model with variable interaction range which is well suited for numerical 
analysis. We decided to consider a ID version of the Random Energy Model (REM) [l^ 



introduced in the first of ref. llj. This consists in a line of mL Ising spins, divided in L 
groups of m spins such that only neighbouring groups of spins interact (thus leading to an 
interaction range of 2m). 

For each group i = 1, we define a state variable taking values 1, ...,2*^. In the 
variables {cTi} the interactions are restricted to nearest neighbors. The Hamiltonian of the 
system is 

L-l 

Hia) = Y,E,{(^^,(^^+l) ■ (1) 

i=0 

For each link the 2^"^ interaction energies Ei{a, r) are quenched random variables extracted 
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from a Gaussian distribution of zero mean and variance 



= m/2 . (2) 

We have considered fixed boundary conditions on the left side (in i = 0), defining ctq = 1, 
and open boundary conditions on the right side (in i — L). In this way we minimize the 
computational effort needed to compute the free-energy Zl, which is expressed as Zl — 
'^^Ziicr), where ZL{a) is given by the recursion relation 

Ze+iia) = J2Mr)e-^'^'^''''\ (3) 

T=l 

with Zq—1 and (3 = 1/T. Computing Zl thus requires 0(L2^'") operations. 

In the m — oo limit, the thermodynamics is simple: the correlations between the en- 
ergy level implied by the one dimensional structure are negligible and, independently of L, 
the free-energy coincides with the one of a REM with 2"*^ states and energies distributed 
according to P{E) oc exp(— £'^/mL): 

^17 /-f-Tlog(2) T>T,, 
hm log Zl^ i (4) 



m— >cx3 



withr,= (2^1^)"'. 



III. THE OBSERVABLES 



Here we define the correlation functions of interest, allowing us to detect a growing static 
length. These are built with the aid of a suitably chosen reference configuration {cr*}i=i,...,L, 
to which one fixes the system outside a window with ^ sites located around the center of 
the system. For convenience we renumber the sites in the central window. Inside 

the window the system is at thermal equilibrium. We investigate the correlation among 
typical in-window configurations B with B* to see whether a characteristic length exist 
such that for window sizes I < ^ — ^* inside the window, while, for I > ^ and a* are 
uncorrelated. 

As detailed in the following, in order to sharpen the transition from correlated to uncor- 
related behavior we decided to fix the reference configuration a* always to the ground state. 
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We then study the thermodynamics of a system which is fixed to the reference configuration 
outside a window of size £: 



cTj = a* Wi < 1 and Wi > i . (5) 

The system has then fixed boundaries and i free variables, {(Ji}i=i^,,,/. Within the window, 
we can define its overlap with respect to the reference configuration as 



q[(T,a ) 

1=1 



7E5K^n- (6) 



Notice that our point-to-set correlation function differs from the one defined in [J] and used 
subsequently which consists in choosing a* as a configuration thermalized at temperature 
T. 

We need some observable estimating the similarity of the typical configuration with re- 
spect to the reference one, and to this end we introduce the following two quantities: 

where the denominator is the "window partition function". The first quantity, Pq, is the 
relative weight of the reference configuration in the window partition function computed at 
inverse temperature /?, while the second quantity, go, is the mean overlap with the reference 
configuration. Both quantities still depend on the quenched disorder and we compute their 



typical values by log(p(yp) = log(po) and log{qtyp) = log(go), where the overline stands for 
the average over the quenched disorder. We expect log(po) and log(go) to be self-averaging, 
since their are related to free-energy differences. 



IV. THEORETICAL ANALYSIS 



In this section we address the problem of an analytic computation of the correlation 
functions. We will first study exactly the asymptotic long range limit m oo. After that 
we will address the problem of finite m effects, that our numerical analysis below reveals to 
be very large. 
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A. The correlation functions for m ^ oo 



The infinite m limit can be understood since in this hmit the correlations between the 
energy level due to the one- dimensional structure of the model become negligible. In this 
case, using this independence approximation, we see that, besides the state a* of energy E*, 
the window has 2™^ — 1 states with energies distributed according to Pe{E) oc 
(there are i sites and i + 1 links!). So that the average density of states is 

Af{E) ~ 2'"^e-^'/™('+^) + 6e,e* (9) 

and the microcanonical entropy per link (divided by m) as a function of the link energy 
e = E/m{i+ 1) is 



-^' + T^M2) H<y/T^log(2) 



Se{e) = { ' ' ' ' V ' (10) 

otherwise 

From this function the canonical thermodynamics can be derived. Before doing that, few 
comments are in order: (a) The constrained entropy is reduced by a constant term with 
respect to the unconstrained case, given by the above formula with £ = oo (see Fig. [1]). (b) 
The choice of the reference configuration a* as the ground state has no effect on the other 
states: the same entropy would be obtained for different choices of a*. Of course the window 
thermodynamics and correlations would depend on the energy of a*. 
If the state a* was absent, the free-energy per link would read 



13 T 



log(2) T > T,Jl + 1/i 



f{P,e) = { 4 i+v^^. • ' (11) 



Including the state a* one therefore has 

/(A £) = min{e*, /(/?,£)}, (12) 



with e* = — A/log(2), that is the ground state energy of the i = oo system. When the 
two terms in Eq. (fT2|) are equal, a first order transition takes place (see Fig. [1]) at inverse 
temperature 

AW^2v/i5i(2)(i-^)-A(i-^). (13) 

which in turns defines a temperature dependent critical length 
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FIG. 1: Microcanonical entropy for the full system = oo, dotted line) and for a window of 
size ^ = 4 (thick line). A system prepared in the ground state of the full system (big dot) and 
constrained in a window of size £ = 4 makes a first order transition at inverse temperature /3c{i = 4) 
(the slope of the thin line). 

separating the confined regime i < ic where ptyp — qtyp — 1 from the deconfined regime 
£> £c where ptyp = qtyp = 0. 

The size of the critical window diverges as expected at the critical temperature, where the 
configurationl entropy vanishes. We find that this critical length is quardatic in the inverse 
of T — Tc] had we chosen the reference state a* with a different rule, the result would have 
been different. For example a direct calculation shows that choosing a* with Boltzmann 
probability at temperature T implies a linear critical length in 1/{T — T^). 

We can understand better the structure of the excitations studying the window free-energy 
as a function of the overlap i.e. the free-energy of configurations that do not coincide with 
B* on exaclty d sites among the £ of the window, with q — 1 — d/i. For simplicity we can 
consider the contribution of "one bubble configurations" where all the d sites in question are 
contiguous. We show below that configurations with more that one bubble are exponentially 
unprobable for large m values. 

For large m the dominant contribution to the free-energy per link f{(3,£,d) — 



7 



B 




lim^ 



FIG. 2: The free-energy barrier to relax from the ground state a*. 
F/m{i+ 1) is independent of the position of the bubble and reads, for d = 1, 



T < T,^l + 1/d . 



(15) 



For d = the free-energy is simply given by /(/?,£, 0) = e*. As one can explicitly see, 
/ is monotonically decreasing in d: the completely open configuration is always the most 
favoured among the ones with d > 1. Notice that at low temperature f{P,i, 1) > f{P,i,0) 
and the difference B{P) = (/ + !) [/(/?, i, 1) - f{P, i, 0)] can be interpreted relaxation 
free-energy barrier for a system prepared in the ground state a*. The barrier B{j3), plotted 
in Fig. [21 is ^-independent and vanishes at a temperature T = T^/ (1 — 1/^/2). 

If we remove the assumption of considering only one-bubble configurations, the free- 
energy in Eg. (1151) becomes 



f{P,i,d) 



min f{f3J,d,b) 

6g{l,d/2} 



(16) 



where f{f3,i,d,b) is the free-energy of configurations differing in d variables from a* and 
having b bubbles, given by the following expression in the large m limit: 



f{P,i,d,b) 



e+l-d-b 



T 



1+h/d 



log(2) T>T,^l + b/d 



aog(2) 



1 



d+b 
£+1 V l+b/d 



log(2) 



T<T,^lTb/d 



(17) 
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It is easy to verify that the minimum in Eg. (1160 is always achieved in 6 = 1, i.e. on one- 
bubble configurations. Multi-bubble configurations can only modify the corrections to the 
leading behavior in m. 



B. Analysis of the Ground State 

We would like to present here some attempts to take into account finite m contributions. 
Corrections to the asymptotic result have two sources: the correlations between the levels 
and sample-to-sample fiuctuations. Though we were not able to deal with the former, we 
could analyze some of the latter. 

Actually we derive some analytical results under 2 main approximations, namely (i) 
energy levels are basically treated as uncorrelated and (ii) the energy of the reference config- 
uration (the ground state energy E*) is considered to be evenly distributed among the links, 
each one having a local energy me* = E*/L (please note that the entire system is made of L 
links, while the window had i + 1 links). We will see below that numerical evidence shows 
that this is the case not too close to the boundaries i = and i = L. 

The distribution of e* is known for L = 1, since in that case me* corresponds to the 
minimum among 2™ independent random Gaussian variables of variance m/2, that is 

^ 1) 4 - VtoiM + + '°«'^^^)' + ^■'^ + ('^] . (18) 

4mv/log(2) \ J 

where X is a Gumbel distributed variable, i.e. F{X > x) = e'^"" . Similarly a closed formula 
can be obtained for L = 2 which corresponds to a two level GREM. Unfortunately as 
soon as L > 2 there are no exact results on the ground state energy of the model. In 
this case an analitical upper bound can be simply constructed by the following greedy 
algorithm: given that ao is fixed, assign (Xi to the value minimizing E{cxo,ai) and repeat 
the procedure recursively on the next variable; at each step the link energy has the same 
probability distribution as me*(L = 1), and so the global ground state energy satisfies 
E*{m,L) < mLe*{L = 1). 

Our numerical data suggest this bound to be tight at the leading order in m for any value 
of L. More precisely, we find numerically that the mean ground state energy can be very 
well fitted, for large values of L, by the following formula: 



E*(m,L) r — — A B 



+ ITTTTT , (19) 
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FIG. 3: The mean value of the ground state hnk energy as a function of the position in the system 
for L = 10 and L = 20. On the left boundary the configuration is fixed, while on the right end the 
system is free. The horizontal line is — ylog(2). 

with A ~ 0.7 and B ~ 0.4. This behavior clearly shows that the convergence to the 
asymptotic intensive energy, — -^log(2), becomes faster increasing L: for L = 1 corrections 
are ( ) and they become 0(m~^/^) in the L — oo limit. We see from Fig.[3]that already 
for L ~ 10, not too close to the boundaries, ground states link energies are independent of 
L, and their numerical values are well represented by the previous formula with L = oo. 
Obviuosly, given the values of m we can study, formula ( fT9l) has to be taken as an empirical 
interpolating function. We find from our data that in ground state configurations, link 
energies have very small sample-to-sample fluctuations, which decrease for larger m values: 
for this reason considering only mean values for the link energies is a good approximation. 
In order to minimize finite L effects and have a homogeneous ground state inside the window 
we find that it was enough to consider sistem sizes L = i + 20, i.e. 10 sites between the 
window and system boundaries. 

Once understood the ground state structure, let us now turn to the estimate of the 
window correlation functions. 
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C. Finite m estimates of the correlation functions 



Under the assumptions stated above, the weight, in the window partition function, of all 
the configurations differing in d variables with respect to a* is given by 

dz e"^"^"^^^!^ 

Z,= (i+l-d) (2™ - l)-^ _ ^-fS{i-d)me' ^ (20) 

/(T+i)m.* dz 

where the first term gives the number of ways to place a bubble of size d in a window of 
size i, the second term counts the number of configurations of the d variables which have 
to differ from a*, the fraction is the average of e~^^ over the p.d.f. of the energies of the 
bubble (it is the sum of d + 1 Gaussian variables of variance m/2, bounded from below by 
the ground state energy, {d + l)me*) and the last term is given by the / — d links having the 
ground state energy. In equation (1201) we have that d G {!,...,£}, while the weight of the 
ground state is given by Zq = exp{—f3E*) = exp{—(3{i + l)me*). We do not write explicitly 
the dependence of Zd on (3, i, m and e* in order to keep the notation light. 

Z^ is an annealed approximation for the window partition function at a fixed distance 
from the ground state. Still, the fact that we keep the dependence on e* explicit is important 
in order to control some fluctuations: e.g. both log(po) and log(go) are given by free-energy 
differences, where the dependence on e* is partially canceled out, and their average over 
e* can be done without any approximation. The two observables we are interested in are 
indeed given by 



Po = ^^, go = , (21) 

and can be easily computed by evaluating numerically the integrals in the definition of Z^, 
once the p.d.f. of e* is known. We have measured numerically such a distribution, but once 
we plugged it into Eq.(l20l) we discovered that the observables we are interested in {ptyp and 
qtyp) mainly depend on the mean of e*, being such a distribution very narrow. Moreover we 
are mostly interested in the dependence of these observables on m in order to understand 
the approach to the m —>■ oo limit, and the average of e* carries the largest dependence on 
m. 

For these reasons the analytical curves we are going to compare with numerical data in 
the next section have been obtained using a non- fluctuating value for e*, give by Eq.( JT9ll 



that is e* = — ^log(2) + 0.726/m'^/^. As we show below, this dependence on the interaction 



11 



range is already enough to produce strong finite m effects. Remind that, in the m oo 
Unait, the logarithm of Zd is given by the free-energy in Eq. ffTSj) . and both ptyp and qtyp 
should drop from 1 to when i crosses the value of idP) given by Eq. (IT^ . 

V. NUMERICAL RESULTS 

The aim of this section is to compute numerically the above defined critical length scale 
for the ID random energy model. The numerical experiment we have performed consists in: 

1. computing the ground state of a system of size L; 

2. fixing the ground state configuration outside a window of size i; 

3. computing qtyp and ptyp in order to see whether there is a first order transition in these 
quantities varying the window size i. 

In the m oo limit we expect such a transition when the window size crosses the value 
ic{f3) given in Eq.f ll4l) . For finite values of m the system cannot have any transition (it 
is one-dimensional), but still the crossover may be very sharp. Our main interest is in 
understanding how much the behavior of finite m systems resembles the mean-field (i.e. 
m —>■ oo) limit and how fast is the convergence. 

As explained in the previous section, we take the size of the system L larger than the 
size of the window in order to avoid boundary effects; that is, to all practical purposes we 
are working in the L ^ cxd limit. 

Thanks to the one- dimensional topology all the experiments can be done exactly by 
tranfer matrix methods. Unfortunately for each link we have a different random matrix 
with 2™ X 2™ entries; for this reason we are forced to small values of m (actually we use 
m = 6,8, 10, 12). Please note that these m values are not so small: the number of degrees 
of freedom per region (2™) is comparable or even larger than the number of particles within 
a typical region studied in reahstic models of glassy systems Since we are interested in 
computing the free-energy at a given value of the overlap with the reference configuration, 
the transfer matrix computation is slightly more complicated and requires a total time of 
order (9(£^2^'"). The average over the disorder is done with at least 1000 samples for any m 
value. 

We are going to present results for temperature T = 0.8, which is a very reasonable value 
(in the m ^ oo limit the critical temperature is Tc = 0.60056...), since the critical window 
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FIG. 4: The overlap qtyp as a function of ^ at T = 0.8. Left panel: numerical values obtained 
throught the transfer matrix algorithm for various values of m. Right panel: analytic curves 
obtained through the approximations discussed in the text. The vertical lines marks the value of 
ic in the m — > cxd limit. 
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FIG. 5: Same as figured] for the quantity ptyp- 

size is 4(T = 0.8) = 15.09. 

In Fig. m and Owe show respectiely log(gjyp) and \og{ptyp) as a a function of i. Left panels 
reports data from exact numerical computations, while right panels show the outcome of 
the analytical approximated analysis. The vertical line is the critical window size ic- 

Some comments are in order. The behavior of all the curves for different m hints at a 
first order transition for m — oo separating a high overlap region at small i from a zero 
overlap region at large i. This behavior is in agreement with the prediction of the mosaic 
theory, however, the convergence is very slow! Indeed so slow, that it does not allow an 
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estimate of the speed of convergence. The crossing point of numerical data for log{qtyp) 
is around i = 10, well below the predicted ic = 15.09. In principle one could argue that 
the one-dimensional model may have a first order transition at a lower value of ic, but the 
crossing point of numerical data for \og{ptyp), taking place around i = 17, suggests that the 
crossing point is strongly dependent on m and converges for m ^ oo somewhere between 
10 and 17 (we are assuming that both qtyp and ptyp have a jump at the same value of i for 
m —>■ oo). 

Still more evident indications of strong finite m effects come from the analytical curves 
(see right panels of Fig. [Hand Fig. [5]): these have been computed from Eq.( l20i) and Eq.( l2T|l 
with e* = -^\og{2) + 0.726 m see Eq. (|T9l) . which is the best interpolation for the 
ground state energy in the window, far from the boundaries. Although these curves have 
been obtained under some approximations, they look qualitatively very similar to the exact 
numerical data, and also quantitatively are not far from the data. For the analytical curves 
we know that they have a jump in i = ic in the m oo limit, still for the present values of 
m they show a crossing point quite far from ic- 

Moreover the value of the overlap at the crossing point may be very small, depending 
on the overlap one is looking at (see e.g. the value of ptyp at the crossing point). For this 
reason may be very difficult to locate the crossing point (remember that our model has a 
very strong random first order transition in the m —>■ oo limit, and most probably things 
work even worst in more realistic models!). 

We remark that simulating the model for a single value of m it would be difficult to 
claim any agreement with the IRSB theory of glasses and the mosaic state: it is comparing 
different values on the interaction range m that the agreement becomes apparent. One could 
argument that one dimension is the worst possibility to observe any behavior reminiscent 
of a phase transition and in higher dimension the situation could be more favourable to 
the theory. Recent simulations of more realistic binary Lennard- Jones mixtures 7[| however, 
failed to identify a sharp mosaic length. 

In order to understand better why qtyp has such strong finite m corrections and show an 
effective crossing point at window sizes smaller than ic, we have studied the window free 
energy as a function of the overlap with respect to the ground state. We show in Fig. [6] such 
a free-energy for T = 0.8, i = 10 and many values of m in order to study the dependence 
on m. We see that, increasing the value of m, all the curves /(g) tends to decrease, but 
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FIG. 6: Window free-energy for £ = 10 and T = 0.8 as a function of the overlap with respect to 
the ground state configuration. 
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FIG. 7: Window free-energy for T = 0.8, m = 12 and many different values of from 6 (top) to 
20 (bottom). 
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corrections to the m — > oo limit are clearly larger for f{q = 1) than for the rest of the 
curve. Please note that f{q = 1) corresponds to the ground state energy, that converges 
in the m — oo limit to — A/log(2) (represented by the horizontal line in the plot). The 
different convergence rate for different q values can be understood also from the analytical 
computation in the previous Section; indeed in the expression for Zd, see Eq. (l20l) . larger 
corrections are for small d values (corresponding to larger q). 

A discrepancy with respect to the analytical computation, is that the free-energy presents 
a minimum at a positive value of the overlap, while in the m ^ oo limit we expect the 
minimum to be in g = 0. This may be one more effect of the slow convergence to the 
mean-field limit. 

In the tentative of extrapolating the numerical results to the m — oo limit, we have fitted 



f{q) data at fixed g, finding that the limit of f{q = 1) is always compatible with — A/log(2), 
while for g < 1 the asymptotic value of /(g) is quite close to that computed numerically 
with m = 12, especially close to the minimum of /(g). 

In Fig.[7]we show the free-energy /(g) for m = 12 (which is very close to the m ^ oo value 
in the low g region) for many i values, ranging from 6 to 20 (top to bottom). The apparent 
first order phase transition between the g = 1 and the small overlap regimes is taking place 
between i = 9 and i = 10 when the minimum goes below f{q = 1) 13(], consistently to what 
we observe in the left panel of Fig. HI Nonetheless, in the m ^ oo we expect the transition 
to take place when the minimum goes below the value — A/log(2), and we see from Fig. [7] 
that this happens around i = 14, much closer to the predicted i^T = 0.8) = 15.09. 



VI. CONCLUSIONS 



The scope of this paper is to study the properties of convergence to the mosaic picture in 
models with larger and larger interaction ranges. We showed that, as it should be expected, 
the behavior of point to set correlations approach the behavior predicted by the mosaic 
picture for large interaction range. The numerical evidence in favor of that comes from a 
differential analysis comparing the behavior for different values of the interaction range m. 
Curves at single values of m do not allow to distinguish mosaic behavior from a single state 
picture where the point to set correlation exhibit a smooth behavior as a function of i. This 
is unfortunate as it indicates that it could be difficult to find confirmations or disprovals of 
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the mosaic picture in realistic glass former models on the basis of the behavior of point to 
set correlations. 

Some papers have recently addressed the study of point-to-set functions in non disordered 
models. Amazingly, the model where the mosaic predictions seems to fit better the data is 
a kinetically constrained model considered in 6|] where a step like behavior of the overlap as 
a function of the window size is observed. Conversely, for a Lennard- Jones binary mixture, 
though it is observed a characteristic length growing with temperature, no step behavior is 
seen. We remark on this purpose that in our data it would been difficult to decide in favor of 
the mosaic picture on the basis of a single value of m. It is only comparing different values 
of m that evidence for the first order jump has been obtained. 

Moreover the convergence to the large m limit is rather slow: finite m curves are very 
smooth and show no precursor of the asymptotic step-like behavior. Ref. [1^ studied a sim- 
ilar Id model with finite interaction range, namely a XORSAT model. The main difference 
with respect to our study is that the model studied in ref. |l2| possesses zero-energy ground 
states and it has been studied only at zero temperature. Despite these differences also in 
ref. 12| large finite-range effects have been found. 

The main effect that we have seen in the model studied here is that a rather sharp 
transition takes place at a finite temperature /5c(£) between a single low-energy ground state 
(i.e. of zero complexity) and a set of higher free-energy states (with positive complexity) 
in a way more or less similar to the mosaic picture. However this transition is plagued by 
large fiuctuations mainly due to the energy of the ground state, which plays a fundamental 
role in determining the critical temperature: the final effect being a sizable smoothing of 
the random first order phase transition at finite value of the interaction range m. 

The conclusions reached in this work suggest that the direct observation of the phase tran- 
sition predicted within the mosaic theory may be rather difficult in realistic models, where 
the interaction range cannot be made very large. A smarter approach for the identification 
of such a transition is likely needed. 
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